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A high temperature expansion is employed to map some complex anisotropic nonhermitian three and 
four dimensional Ising models with algebraic long range interactions into a solvable two dimensional 
variant. We also address the dimensional reductions for anisotropic two dimensional XY and other 
models. For the latter and related systems it is possible to have an effective reduction in the 
dimension without the need of compactifying some dimensions. Some solutions are presented. This 
framework further allows for some very simple general observations. It will be seen that the absence 
of a "phase interference" effect plays an important role in high dimensional problems. A very 
forbidding purely algebraic recursive series solution to the three dimensional nearest neighbor Ising 
model will be given. In the aftermath, the full-blown three dimensional nearest neighbor Ising model 
is exactly mapped onto a single spin 1/2 particle with nontrivial dynamics. All this allows for a 
formal high dimensional Bosonization. 



I. INTRODUCTION 

In this article, a high temperature expansion is em- 
ployed to examine mappings between problems of vari- 
ous dimensionality. By these reductions we will exactly 
solve an anisotropic two dimensional XY and some other 
models. It will be seen that the absence of a "phase in- 
terference" effect plays an important role in high dimen- 
sional problems; this observation will allow us to examine 
various high dimensional models by averaging over newly 
constructed one dimensional systems. In turn, this aver- 
aging process will formally allow us to effectively bosonize 
high dimensional theories. 

In section ||, I am going to start off by discussing a 
momentum space high temperature expansion. Next, by 
analytically continuing the momenta and making them 
complex (section III) , quaternion (section IV), and of a 
general matrix form we will be able to exactly solve a va- 
riety of two, three, and four dimensional Ising and 0{n) 
models. Many of these models are artificial and corre- 
spond to non hermitian anisotropic long range interac- 
tions. Some of the resulting models will be short ranged 
and of real hermitian form. We will exactly solve some 
special simple two dimensional anisotropic XY models 

We will then extend these techniques to the quantum 
case and discuss "2-1-1 dimensional Bethe Ansatz" solu- 
tions in section(^ where we will also solve a special two 
dimensional anisotropic XY model. 



After solving these models, in subsection (VIA) we will 
apply the high temperature expansion to reduce the three 
dimensional nearest neighbor Ising ferromagnet into an 
illusory formal solvable expression written in a determi- 
nant notation. In this method, the high temperature 
coefficient to each order is arrived at by looking a single 
linear relation originating from a matrix of (super expo- 
nentially) increasing order. 

In subsection( VI B), we will map the nearest neighbor 



d dimensional Ising models to an ensemble average over 
a spin chain with long ranged "Coulomb" like interac- 
tion. In this scenario each spin is effectively composed 
of d constituent charges each generating its own shifted 
"Coulomb" like field. This in turn will enable us to map 
the d dimensional Ising model onto a a single spin 1/2 
quantum particle in 0+1 dimensions. We will apply the 
dimensional reduction techniques to quantum spin prob- 
lems where an effective one dimensional fermion problem 
will result. 



In section( VII), I will discuss a "permutational" sym- 
metry present in the 0{n —^ oo) models but violated to 
0{f3*) in finite n problem. This will again allow us to 
map any large n problem onto a one dimensional system 
but this time without the need of averaging. 

One of the main motifs of this article is the absence or 
presence of "momentum interference" effects in high and 
low dimensional problems respectively. 



II. A MOMENTUM SPACE HIGH 
TEMPERATURE EXPANSION 

We will for the most part examine simple classical spin 
models of the type 

1 v^. 



2 ^ 



Vix,y)[S{x)-Sm. 



(1) 



Here, the sites x and y lie on a (generally hypercubic) 
lattice of size N. The spins {S{x)} are normalized and 
have n components. The kernel V{x,y) is translation 
invariant. The n — 1 case simply corresponds to the 
Ising model. 

Later on, when we will discuss quaternion models we 
will replace the scalar product by a slightly more com- 
plicated product. 



1 



The Hamiltonian in the (non-symmetrical) Fourier 
basis is diagonal (/(fc) = ^(^)e"*'''^; F{x) = 



jr Efe /(fc)e*'=-^) and reads 



(2) 



where v{k) and S{k) are the Fourier transforms of V{x) 
and S{x) respectively. 

For simplicity, we will set the lattice constant to unity- 
i.e. on a hypercubic lattice (of side L) with periodic 
boundary conditions (which we will assume throughout) 
the wave-vector components ki = ^jp- where r; is an 
integer (and the real space coordinates xi are integers). 

In the up and coming V = (3V . For an invertible 
V{x^y)^ and < < oo (via a Hubbard-Stratonovich 
transformation [Q): 

Z = TreM-\Y.^i.S.y)S{S)S(,y)} 



= ydct(-A) j ^d^{x)eM~H{r^{m 



(3) 



where 



(4) 



— ln[cosh 77(2?)] 

X 

and analogously for the 0{n > 1) model 



■^ln[ r(n/2) ij^r/'-' In/2-i{\m\) ] (5) 



with /ri/2-i(-z) a Bessel function. 

The cosh?7 or Bessel terms may be viewed as terms 
of constraint securing the normalization of S{x) at every 
site X. 

If a magnetic field h{x) were applied, the argument 
of the cosh would be replaced cosh(r7 -I- h) with a simi- 
lar occurrence for the Bessel function appearing for the 
0{n > 1) models. 

The g— state Potts model can be viewed as a spin model 
in which the q possible polarizations of the spins lie at 
the vertices of a (g — 1) dimensional tetrahedron- in this 
manner the scalar product amongst any two non iden- 
tical spins is cos^^(— l/(g — 1)). Employing this repre- 
sentation, the sum Ej In cosh 77 (if) may be replaced by 
Ex -^'^Ei=i 6xp[S^i • ff\ where the sum is over the q polar- 
izations of the {q — 1) dimensional spin Si. 



The partition function of the Ising spins reads 



Z=^dct(-A)|i?,e- 

X 

^r]{x)V^^{x,y)r]{y) 



JJcosh?7(x). (6) 

s 

For a translationally invariant interaction V{x, y) = 
V{x — y), this may be written in momentum space 



1 

2N 



HE 

X m—0 



772™ (x) 

(2m)! ' 



(7) 



where v{k) is the Fourier transform 



z;(A?) = ^F(x-y)e*'^-^ = /3v(fc). (8) 

X 

where the sum is over all lattice sites x. Thus the parti- 
tion function is , trivially, 

mi,...mjv i—1 ^ ' 

where ( )o denotes an average with respect to the un- 
perturbed Gaussian weight 



exp 



^Y^^-\k)^{k)^{-k) 

k 



(10) 



and A/" is a normalization constant. Each contraction 
{il{x)r]{y)) (or {■q{k)ri{—k))) leads to a factor of V{x, y) — 
l3V{x, y) (or to Pv{k)/N in momentum space). Thus the 
resultant series is an expansion in the inverse tempera- 
ture (3. In the up and coming we will focus attention on 
the momentum space formulation of this series. All the 
momentum space algebra presented above is only a slight 
modification to the well known high temperature ex- 
pansions usually generated by the Hubbard Stratonovich 
transformation directly applied to the real space repre- 
sentation of the fields r]{x) ^j^. It is due to the naivete 
of the author that such a redundant momentum space 
formulation was rederived in the first place. However, as 
we will see, in momentum space some properties of the 
series become much more transparent. To make V in- 
vertible and the expansion convergent, we will shift v{k) 
by a constant 
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v{k) -> v{k) + A 



(11) where in the last sum m'^^„ ^ ml 



such that v{k) is strictly positive. In real space such a 
constant shift amounts to a trivial shift in the on site 
interaction (or chemical potential) 



Vix.y) V{x,y) + A.Ss.y 
For asymptotically large \r]\ the nontrivial 

i/f""' =^ln|coshr?(f)| 



(12) 



(13) 



is linear in jryj. The Gaussian generating 

= ^ E vi£)V-Hx, y)v{y) (14) 



is a positive quadratic form and dominates at large rj. 
So far all of this has been very general. Now let us con- 
sider the d dimensional nearest neighbor Ising model. To 
avoid carrying minus signs around, we will consider an 
Ising antiferromagnet whose partition function is, clas- 
sically, identically equal to that of a ferromagnet. The 
exchange constant J will be set to unity. The correspond- 
ing momentum space kernel 



vik) = A + J2 



(15) 



1=1 



As the reader might guess, the uniform shift A only 
leads to a trivial change in Z (i.e. to a multiplication 
by the constant exp[—(3AN]). In the forthcoming we will 
set A^O. 

Owing to the relation 



dk e'''^ = 2tt5n 



(16) 



all loop integrals are trivial. A simple integral is of the 
type 



M 



iM 



(17) 



(when n = 1 this integral is 2d etc). A related integral 
reads 



d'^k 
{2tj:Y 



n 



i<M 



etc. 



Note that the high temperature expansion becomes sim- 
pler if instead of symmetrizing the interaction (each bond 
being counted twice- once by each of the two interacting 
spins) one considers 



(18) 



){k) — 2y~^ exp[ifc;]. 



(19) 



Here each spin interacts with only d (and not 2d) of 
its neighbors separated from it by one positive distance 
along the d axes ef^^. The factor of 2 originates as each 
bond is now counted only once and therefore the corre- 
sponding bond strength is doubled. The partition func- 
tion Z is trivially unchanged. 

In the most general diagram the propagator momenta 
{<la\ (1 < a < number of propogators) are linear combi- 
nations of the independent loop momenta {fcb};,=*i^^ 

qa = Mabkb- (20) 

Symmetry factors aside, the value of a given diagram 
reads 



d-'kb 



n n(2/3)Sf=i exp[zS,A4,fc^] 



_ ^2/3)P''oP^g^tors 



loops 



n (/4^^trexp[^E^^-^^ 

6=1 y I a,b 



(21) 



- i.e. upon expansion of the outer sum, just a product 
of individual loop integrals of the type encountered in 
Eqn.(|l|) for each 

The latter integral in d = 1 reads 



loops 

n 

6=1 



^(^propagators^^^^ 



(22) 



The situation in d > 1 becomes far richer. One may parti- 
tion the integral into all possible subproducts of momenta 
corresponding to different spatial components I satisfy- 
ing the delta function constraints. In a given diagram 
not all propagators need to correspond to the same spa- 
tial component fc'. The combinatorics of "in how many 
ways loops may be chosen to satisfy the latter delta func- 
tion constraints?" boils down to the standard counting 
of closed loops in real space. 

Some simple technical details of lattice perturbation 
theory vis a vis the standard continuum theories are pre- 
sented in iiotl. 
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III. AN EXACT SOLUTION TO A THREE 
DIMENSIONAL ISING MODEL 

Theorem: 

If the spin spin interaction on a cubic lattice has a non 
hermitian kernel of the form 



1 sinhTrA 



HSx2,l + 6x2-1) 



+ix3{S, 



2:2,1 



'>X2-1) 



(23) 



then, for all real A, the Helmholtz free energy per spin 
a{h = 0,T) is 



Pf = - ln(2 cosh 2/3) 
— j (i0 In - ( 1 + \J 1 — K-^ sin 



(24) 



where 



(25) 



_ 2sinh2/3 
^ ~ cosh^ 2/3 

and the internal energy per spin 

u = -coth2/3[l + -K'i^i(K)j (26) 
where Ki{k) is a complete integral of the first kind 



and 



k' = 2tanh2(2/3) - 1. 



(27) 



(28) 



The proof is quite simple. The expressions presented 
are the corresponding intensive quantities for the two di- 
mensional nearest neighbor Ising model . 

Let us start by writing down the diagrammatic expan- 
sion for the two dimensional Ising model and set 



«(?!, 92) v{qi,q2 + iXqs) 
in all integrals. This effects 



dh 



Jk2N'''°" 



'26 



dk 



26 



and thus an additional factor of L is introduced for each 
independent connected diagram. Thus the free energy 
per unit volume for a system with the momentum space 
kernel v{ki,k2 + iXk^) is the same as for the two di- 
mensional nearest neighbor ferromagnet with the kernel 
v{ki, ^2)- Fourier transforming v{ki, k2 -l-iAfca) and sym- 
metrizing [V{x) + V{—x)]/2 V{x) one finds the com- 
plex real space kernel V{xi, X2, X3) presented above pl| ] 
Q.E.D. 



Corollary 

For the more general 



V{xi,X2,X3) = 

we may define 



-'6x2fi6x3,0{6xiS + Sxi,-l) 
{{6x2S + Sx2-l) 



A2 



+ i{Sx2-l - Sx2.a)-J-} 



ttJ 



2A sinhTrA 



(32) 



(33) 



to write the free energy as the free energy of an 
anisotropic two dimensional nearest neighbor ferromag- 
net with the kernel 



V{xi,X2) = ~Sx2,oi^xi,l + '5x1-1 

-r6xi,oiSx2,i + Sx2-l) 



(34) 



i.e. a nearest neighbor Ising ferromagnet having the ratio 
of the exchange constants along the xi and X2 axes equal 
to r = J2 / Ji (throughout we set Ji to unity) . Introduc- 
ing 



1 



sinh 2/3 sinh 2/3r 



(35) 



and defining 



F{e) = ln{2 [cosh 2/3 cosh 2/3r 
+u-\l - 2ucose + u^y/^]}, 



(29) 

we may write the free energy density as 



dfeae-^'^^^""" (30) 



(36) 



(37) 



for all individual loop integrals. The integrals are nonva- 
nishing only if {/Vb = 0} for all loops b. When iV'°°P = 
for a given loop the second (fca integration) leads to a 
trivial multiplicative constant (by one). In the original 
summation 



Note that the same result applies for the "ferromag- 
netic" 



fet 



(27r)'* 



V{xi,X2,X3) = 
— J^xifi 



d'^kh 



(31) 



"'5x2, 0*5x3, 0('5xi4 + '5x1, -1) 



-'i('5x2 



'5x 

,1 + '5x2, -1) 

'5^2,i)yl- 



(38) 
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The proof is simple- the partition function is identically 
the same if evaluated with the kernel in Eqn.(|3^) instead 
of that in Eqn.(p2|) if the spins are flipped 

(39) 



Six) ^ {-irS{x). 



What about correlation functions? The equivalence of 
the partition functions Z[f3, h{xi, X2)\ implies that in the 
planar directions the correlations are the same as for two 
dimensional ferromagnet. 

Changing the single momentum coordinate 



effects 



v{ki,k2) v{ki, k2 + iXks) 



G{ki,k2) G{ki,k2 + iXk3 



(40) 



in the momentum space correlation function. The alge- 
bra is identically the same. 

The connected spatial correlations (for T ^ Tc) are 
exponentially damped along the X2 axis as they are in 
the usual two dimensional nearest neighbor ferromagnet 
(see for the two dimensional correlation function). 

Of course, all of this can also be extended to one- 
dimensional like spin systems. 

The two dimensional spin-spin kernel 



V{xi,X2) = ( 



sinhTrA 

'^1 TT'. 2"a:i,0 



X^ + xi 



^(^2,1 + '^2:2,-1) 

ix2{Sx2S - Sx2-i)] (41) 



leads to one dimensional behavior with oscillatory oscilla- 
tions along the X2- More precisely, the momentum space 
correlator for any 0{n) system in one dimension reads 



G(fc) = 



1 



1 



(42) 



For the Ising (n = 1) system r = — [lntanh/3]. One 
may set k ~ ki + iXk2 and compute the inverse Fourier 
transform. The correlations along X2 are oscillatory with 
temperature dependent wavevectors. In other words, the 
effective correlation length along X2 is infinite. 

If the two dimensional spin-spin interactions in 
Eqn. (|4l|) are augmented by an additional on-site mag- 
netic field h then the two dimensional partition function 
reads 



Z = X^ + X 



>7V 



where, as usual. 



cosh Ph ± \Je^P sinh^ /3/i + e-2/3. 



(43) 



(44) 



That "complexifying" the coordinates should not 
change the physics is intuitively obvious: if one shifts 



fci — > fci -f Xk2 = k'l with real A in the one dimensional 
v{ki) then the resulting kernel is v{k'i) plainly describes 
a stack on chains parallel to (1,A); the spins interact 
along the chain direction yet the chains do not interact 
amongst themselves. The free energy density should be 
identically equal to that of a one dimensional system- no 
possible dependence on A can occur. All that we have 
done in the above is allow A to become complex. 
In 



v{ki,k2) — > v{ki,k2 + Afca) 



(45) 



the correlation functions along a direction orthogonal to 
the direction (1,A) in the (0:2, a^s) plane vanish. If A is 
complex then there are no correlations along an orthog- 
onal direction in the "complex" space. 

Formally, all this stems from the trivial fact that the 
measure dkidk2 = const{dkdk*), whereas v{k) depends 
only on the single complex coordinate k. 

Other trivial generalizations are possible. For example, 
we may make both fci and k2 complex in v{ki,k2) to 
generate a four dimensional spin kernel 



V{XI,X2.X3,X4) 

1 // 1 Nro SinhTrAl „ ^^ / r c- N 

'-((-1) X2 I 2 '^=^1.0 Al(^X24 +'5^2,-l) 

+1x3(6x2^ - ^0:2 -1)] 



^^ sinh^Aa 



A2 -\- 



X2i(>XiA 



+1X4(6x3^1 - Sx3-l)] 



(46) 



which leads to the canonical (nearest neighbor) two 
dimensional behavior for arbitrary Ai and A2. 

All of this need not be restricted to cubic lattice mod- 
els. We may also analytically continue v{k) of the trian- 
gular lattice etc. 

For the triangular antiferromagnet/ferromagnet 



fci Vs 

v{k) = ±2[cosfci -I- cos(— -I- -^^2) 



+ cos(^ 



fci Vs 



k2)]. 



(47) 



Making k2 complex leads to interactions on a layered 
triangular lattice etc. 



IV. HIGHER 0{N > 1) MODELS 

The multi-component 0{n > 1) spin models display 
a far richer variety of possible higher dimensional exten- 
sions. 

A need not be only a complex number, it may also be 
quaternion or correspond to more general set of matrices. 
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Let us examine the extension to "quaternion" k (or 
with the matrices icti, icr2,and ia^ (with {aa} the Pauli 
matrices) taking on the role of i,j and k). If we were 
deahng with an 0(2) (or XY) model, the same high tem- 
perature expansion could be reproduced: 

We may envisage a trivial extension to the usual scalar 
product Hamiltonian 



H 



x-y)S{x) ■ S{y) 



(48) 



to one in which the kernel is no longer diagonal in the 
internal spin indices a, (3 — 1,2. 



H = ^^Va,/3{X - y)Sa{x)S(3{y). 



(49) 



The high temperature expansion may be reproduced 
as before with the kernel Va,/3{x — y) being generated by 
each contraction of ija^x) and rj/sijj), or in momentum 
space Va.p{k)5^j^p ^ is generated by the contraction of 

r]a{k) and rjp^k'). 

Let us start off with a nearest neighbor one dimensional 
XY chain and consider the transformation 

exp[ifci] exp[i(/ci + iAi/c2f i + 1X2^71,(72 + iX^k^a^)] 

= exp[ifci — w ■ a]. (50) 

In the argument of the exponential the identity matrix 1 
commutes with w ■ a and the exponential may be trivially 
expanded as 



exp[ifci " w ■ 3] = exp[«fci] exp[— • 3]. 



(51) 



Once again in higher dimensions (this time d — A) the ki 
integration will reproduce the familiar 



I dki exp[iNki] — 2n5N^o 

J — TT 



(52) 



for each independent loop momentum k and the integra- 
tion over the remaining k2,k3, and k^ components will 
simply yield multiplicative constants. 
Thus the kernel 



v{k) = exp[ifci — w ■ a] + h.c. 



(53) 



will generate a four dimensional XY model which has 
the partition function of a nearest neighbor one dimen- 
sional XY chain. The corresponding real space kernel 
V'(a;i,a;2,X3,a;4) 



Vaf3ixi,X2,X3,X4) 



7T pTT .4 



d'^k 



(27r)4 6^P[*^' ^](6^pW=i 
+iXik2ai + iX2ki,(T2 + iXzkia-i)]aii + h.c^ . 



In general, this kernel is no longer diagonal in the internal 
spin coordinates (unless w happens to be oriented along 

CTs). 



1=2 

Aifc2 A2fc3 Xzki 

= (-r^>-r^,-r^) (54) 
\w\ \w\ \w\ 

and 

exp[— (T • w] — cosh(|i(J|) — sinh(|it;|)((T • w) (55) 

For general {Xi\^^^ the inverse Fourier transform is non- 
trivial. The scalar product is replaced by a more compli- 
cated product amongst the components. The real space 
interaction kernel is once again algebraically long ranged 
along the X2,X2 and x^ axes. 

In this manner we may also generate an infinite number 
of solvable models with a real hermitian kernel. We 
may state another 

Theorem 

The hermitian kernel 



v{k) — —2 cos(fci + k2(T'i) 



(56) 



leads to the partition function of a one dimensional 0(2) 
chain. 

The proof amounts to a reproduction of the calcula- 
tions given above for the four dimensional 0(2) system. 

In real space the latter kernel (Eqn.(^)) leads to the 
0(2) Hamiltonian 

H^- S^sf- E (57) 

{ij) along (1,1) (ij) along (1,-1) 

where the n = 1 (or x component) of the spins inter- 
act in the first term and only the y-components of the 
XY spins appear in the second term; the indices i and j 
are the two dimensional square lattice coordinates. The 
first term corresponds to interactions along the (1,1) di- 
rection in the plane and the second term corresponds to 
interactions along the (1,-1) diagonal. The two spin 



components satisfy [S, 



(l)l2 



(2)l2 



1- i.e. are normal- 



ized at every site i. By our mapping, this model trivially 
has the partition function of a one dimensional nearest 
neighbor 0(2) spin chain. 

Similarly, the partition function for the three dimen- 
sional kernel v{k) — — 2cos(fci -1-^2-1- k-^a^) is that of 
a two dimensional nearest neighbor XY model exhibit- 
ing Kosterlitz-Thouless like behavior. Here each of two 
spin components interacts within a different subplane. 
The coupling between the two spin components due to 
the normalization constraint apparently plays no role in 
leaving the system two dimensional. 
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The reader can easily see how a d! dimensional system 
can be made have an effective dimensionality 1 < d < d! 
without the physical need of actual compactification. 

Polarization dependent hopping (or "interactions" ) 
akin to those in Eqn.(|57|) may be of relevance in dis- 
cussing directed orbital problems in two dimensions. 



V. SPECIAL TWO DIMENSIONAL BETHE 
ANSATZ SOLUTIONS 

Let us also write the general Hubbard Stratonovich 
transformation for a Quantum system 

Consider a system of fermions with the Hamiltonian 



(58) 



Following Negele and Orland [g| it is convenient to in- 
troduce 



KaS = TaS — - ^ Va 



L 

2 



1/3/3(5 



(59) 



/3 



where, for our applications, a, f3, 7, and S are two dimen- 
sional position vectors. 

For this time independent Hamiltonian 



exp[-iH{tf - U)] 

TT i f*^ 

[[DT]ap{t)exp[- / dtr]af3it)vaf3'rSV'rsit)] 

Q/3 



(60) 



Again a momentum space expansion is possible in 
Vapit)- If an exact Bethe Ansatz solution is known for 
nearest neighbor hopping then in some cases it may be 
extended to two dimensions where the hopping matrix 
element K{xi,X2) for separation x = (3:1,2:2) — a — P 
attains exactly the same form as the two-dimensional 
kernel V{xi,X2) just given previously. More generally, 
the Trotter formula may be employed. 

Loosely speaking, if in a fictitious electronic system the 
hopping matrix element would be the non hermitian 

rx^ \ 1/^/ iNa;, sinhTrAi 
Kix„X2) = --[i-l)^^^6.,^o 

^1(^x2,1 + 1^x2,-1) 

+ 1X3(6x2,1 - Sx2-l) ^ (61) 



then the system would essentially a one dimensional Lut- 
tinger liquid along xi. If the kernel would instead decay 
algebraically along the imaginary time axis the problem 
would become that of dissipating system. 

A hermitian hopping matrix element could also do the 
trick if we were to consider polarization dependent hop- 
ping (or "interactions") analogous to those in Eqn.([57|) 
and slightly more complicated variants. As noted earlier, 
such models may be of relevance in discussing directed 
orbital problems in two dimensions. 

For mapping the two dimensional S — 1/2 Quantum 
problem of Eqn. (p7|) onto a nearest neighbor XY chain 
we find that the free energy density is given by 



/?/ = - / dx ln[l + exp(4/3cos(27rx))]. 



(62) 



VI. PHASE INTERFERENCE AND A MAPPING 
ONTO A SINGLE SPIN PROBLEM 

We will now map high dimensional (d > 1) problems 
onto translationally invariant problems in one dimension. 

For the nearest neighbor ferromagnet, the basic idea is 
to make a comparison between the the kernels 



Vd — — 2y~^cos k, 
1=1 

d 

vi = — 2^cosq/c] 



(63) 



1=1 



in d and in one spatial dimensions respectively. Note that 
the one dimensional problem has, in general, d different 
harmonics of the single momentum coordinate fci . If the 
coefficients c/ are integers equal to Ri then the one di- 
mensional problem amounts to a ferromagnetic chain in 
which each spin interacts with d other spins at distances 
{Ri}f=i away. When evaluating the loop integrals, we 
will find that the partition function/ free energy deviate 
from their d dimensional values due to "interference" be- 
tween the various harmonics in ■Oi(A:i). If all harmonics 
acted independently then the d dimensional result would 
be reproduced. The advantages of the form vi are ob- 
vious. Perhaps the most promising alley is that of "in- 
commensurate dimensional reduction" . By this method, 
we may examine the d-dimensional problem described by 
Vd by a bosonization of the one-dimensional problem of 
■Oi. A solution to the one dimensional problem posed by 
■Oi followed by an average over the coefficients {c;} will 
immediately lead to the corresponding d— dimensional 
quantities. 
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A. Commensurate dimensional reduction 

Consider a one dimensional ferromagnetic lattice 
model with a nearest neighbor (or Range — 1) interac- 
tion and a Range = n interaction: v{ki) = — 2(cosfci + 
cosn/ci). The loop integrals containing terms of the 
cos fci origin and terms stemming from cos nki are in- 
dependent to low orders. Mixed terms can survive only 
to orders (l/r)"+^ and higher. The delta function con- 
straints '^qi ~ Q would generate terms identical to 
those of the two dimensional nearest neighbor model with 
^(fci, — — 2(cosfci -I- cosfc2) with the two delta func- 
tion constraints for the two separate components of the 
momentum, '^qi = ^ and 92 = 0, at every vertex. 
Thus, in a sense, this system is two dimensional up to 
T^". Similarly, if a one dimensional chain has interac- 
tions of length one lattice constant, R = n and R = m 
i.e. 

v{ki) — — 2(cosfci -t- cosnfci + cosmfci) (64) 

then this system, as fleshed out by a 1/T expansion for 
the partition function or for the free energy, is three di- 
mensional up to order ~ n \i m ~ r? . (If m = n + 1 
then the triangular ferromagnet is generated). Similar, 
higher dimensional extensions similarly follow from the 
lack of commensurability of the cosine arguments (for 
four dimensions 

v(k\) = — 2(cosfci + cosnfci -|- cosmfci -I- cossfci) (65) 

with s = n^ etc.) 

The corresponding transfer matrices are trivial to write 
down. For the three dimensional case: 

(si...s„2|T|si...s^j2) = exp[/3{i(siS2 + ... + s„2_is„2) 

+ ^K4 + - + <2_i<2) 

+ -{SiSn + l + S2Sn+2 + ■■■ + S„2_„.S„2) 

+ 2(*i*n+l + ^2*'i+2 + ••■ + •''n2_„S^2) 
+ Sn2s'i + (s„2_„^is']^ + S„2_„_|_2S2 + ... + Sn^sJi) 
+ {sis\ + ... + S„2S^2)}] 

= exp[/37;,]. (66) 

In higher dimensional generalizations, Tji will be slightly 
more nested with the span n^ replaced by n"*^^. In 
Eqn.(|6^), the matrix indices i and j are written in bi- 
nary numerals in terms of the n^ spins 

j = 5](s„-Hl)2"-2 (67) 

a=l 

where a is the one dimensional coordinate along the bra 
spin indices. By Eqn.(|66|), Tji are expressed in terms of 



a sum over two digit products where the digits are those 
that appear in the binary (length n spin) representation 
of the coordinates j and i. 

As usual, the partition function 

Z = Tr[T^] = ^ Af (68) 

i 

where {A^} are the eigenvalues of the transfer matrix. In 
the following Amax will denote the largest eigenvalue. 

Note that the transfer matrix eigenvalues correspond 
to periodic boundary conditions as strictly required: we 
employed translational invariance to write the partition 
function expansion in Fourier space. If the boundary 
conditions are not periodic all this is void. 

By the lack of interference effects, the connected dia- 
grams for the free energy are correct to O (/?"■) for inter- 
actions with m = n^ (Eqn.^4|). This along with 

p 

(where the last equality holds in the limit of large sys- 
tem size N) and a canonical power expansion having an 
infinite radius of convergence 

A.na. = exph E fpP'\ - E (ii E fvP'r) 

p m— ' p 

= E&p/5'' (70) 
p 

imply that A^ax is correct to O (/?"). 

If we already know the lower order coefficients corre- 
sponding to 0(/3P) with p = 0, 1, 2, 3, (n — 1) when we 
consider the Range = n problem, then we may set out 
to compute is the coefficient of This will amount to 
extracting only a single unknown (i.e. the term corre- 
sponding to the coefficient of /3") from a transfer matrix 
eigenvalue equation. This can be done to the next order 
etc. recursively. A simple yet very long single linear rela- 
tion gives the largest transfer matrix eigenvalue to each 
higher order in the inverse temperature /3. 

Explicitly, we can write longhand 

^2 
2" 

det(T - A) = e,,...,^^, l[[T,^,{f3) - A(/3)<5,,J = (71) 

All one has to do is to yank out the 0{(3") term from 

2 

the sum over products of 2" terms. If we pull out a 
power pji- from each of the elements [Tji. {(i) — \{P)5ji.] 
then those will need to satisfy 

2 

2" 

E^ji. = (72) 



8 



In each element of the transfer matrix, exp[/97ji], the 
coefHcient of is triviaUy Tj'^/pl. For each permutation 
(or "path" of coordinates within the matrix) to be 
summed for the evaluation of the determinant 



path 



E 



E 



n 



(73) 



is exactly the product of the coefficients of (3^^^ such that 
those lead a net power /3" (i.e. satisfying Eqn.(|7^)) from 
that particular path. 

It follows that net contributions stemming from Tji. (/3) 
in Eqn.(|7l|) are 



^path [ 



W (74) 



all paths 



(jij) in a given path 



where Cpath simply denotes the sign of the given permu- 

tation {ij}. 

To take into account the products including A(/3) let 
us define 



D 



(ti;p) _ 
2m 



1 



{n-py. 



E 



all 2m paths 



^pathl ^ ^ ^ 

(jij) in a given path 



(75) 



where the summation is over paths going through all 
possible {2m) given points on the diagonal and de- 
notes a summation over Tji^ in the given paths sans the 
contributions from the {2m) diagonal points. The de- 
terminant paths in J^aii 2m paths '^^^^ p3.ss through more 
then (2m) diagonal points- it is just that we need to sum 
over all those that in their pass also traverse all possible 
sets of (2m) given points on the diagonal and where those 
points are excluded from the second Tji^ summation. 

The /3" component of Eqn.([7l|) is 

FM+J2J2DtP^ E ^Pi-^p--0 (76) 

where {&p}p^Q appear in the expansion of the largest 
eigenvalue (Eqn.([70|)). 

If the coefficients {bp}plQ then Eqn.(|7^) where &„ ap- 
pears there only linearly gives 6„ and one may then pro- 
ceed to look at the eigenvalue equation for the next value 
of n. 

In the absence of an external field, for n — 1 the 
largest eigenvalue for the d dimensional problem is A_|_ = 
2 cosh ci/3 and other eigenvalue A_ = 2sinhd/3. For 
the three dimensional case setting m = = n = 1 
leads to the one dimensional momentum space kernel 
v{ki) = — 6cosfci which in real space corresponds to 



H 



4+1- 



(77) 



Of the two eigenvalues only A-|_ has a nonzero bo- Know- 
ing the values bo = 2 and fei = of the largest eigenvalue 
we may proceed to find 62 from the Eqn. (^6|) when n ^ 2 
etc. This, of course, may be extended to systems with 
magnetic fields. 



Setting 60 = 2 in Eqn. (|76[) we find 



n-1 

E E bp^^p^^i 

p=0 p=pi+P2 



(n;p) 



+E E 



p=0 m=2 



E 

^^."^^Pi— P with PiT^P 



bpi ■ --^Pn 



(78) 



This is the rather explicit recursive relation giving 6„ 
once {bp}p~Q are known. To find out the coefficient of 5„ 
let us look at 



D 



{n;n) 



^ ^ ^pathi 
exactly 2s diagonal paths 



(79) 



where the summation admits only those permutations 
that lead to exactly 2s diagonal elements. For the 2" 
diagonal path in the matrix T: Spath = 1- For all paths 
threading 2" ^ diagonal sites epath = — 1 and the num- 

ber of such paths is {^^ ). For the general case the last 
quantity amounts to 



D 



n;n) 



2s 



(80) 



Evaluating, we find 



2" 



(2" ) 



s— m 

m 
2" 



^Fp,[{l,^-2-'-\;, + ^^2-'-% 



(4m)!(2" - 4m)! 

1 m 

2 + Y 



(81) 



where Fpq is a generalized hypergeometric function. 

Insofar as simple geometric visualization is concerned, 
it is amusing to note, as shown in Fig.(0) for the two 
dimensional case, that employing the usual high temper- 
ature expansion in powers of tanh /3 one would reach the 
same conclusion regarding the correctness (up to order 
n) of the partition function evaluated for a width n slab 
vis a vis the partition function of the two dimensional 
system. We may look at the a finite thickness {n) slab of 
the two dimensional lattice along which we apply peri- 
odic boundary conditions. Let us now draw a string going 
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along one row of length n, after which it would jump to 
the next row, scan it for n sites, jump to the next one, 
and so on. On the one dimensional laced string, the sys- 
tem is translationally invariant and the interactions are of 
ranges R = 1, n. Or, explicitly, by counting the number 
of closed loops in real space (employing the standard, 
slightly different, diagrammatic expansion in powers of 
tanh (3) , we see that the terms in this expansion are also 
identical up to order tanh" /3 . 

A lacing of a two dimensional slab by a one dimensional 
string, as shown in the figure, is one of the backbones 
of Density Matrix Renormalization Group Theory when 
applied to two dimensional problems. 



B. Incommensurate dimensional reduction 



Theorem: The one dimensional real space kernel 

d 

V{xi) = -2^[(c/ + xi) sin7r(c; - xi) 
1=1 

+ {ci - xi)] sin7r(a;i + c;)] 
x[ici ^ xi)ici + xi)]-\ (82) 

(i.e. a scenario in which each spin is effectively composed 
of d shifted "Coulombic sources" (more precisely, each 
spin is composed of d pairs of "charges" generating sine 
potentials)) will give rise to the exact d dimensional near- 
est neighbor partition function and free energies when 
those quantities are averaged over {ci}f^-^. 

Once again the proof is not too involved. The basic 
idea is that in momentum space this will give rise to 



v{ki) — — 2y~^ cos Cjki 



(83) 



1=1 



which when averaged over incommensurate {c/} will an- 
nihilate all phase coherent terms and reproduce the ex- 
pansion with the d dimensional kernel in Eqn. ([l5|) . 

Many different measures for dfi{ci) can be chosen. Per- 
haps the simplest one is 



(84) 



2A ' 



for each of the d coefficients {c/}. 

Performing the {ci}f^-^ integrations after the loop inte- 
grals over {ka} we will find that the most general integral 
is of the form 



1=1 



dkb 

b a 1=1 



A one dimensional 



parameterization of the (n x L) slab 



3n - 
• 


f 1 

• 


mm, m 


An - 1 
m , 


271- 


n 






n + 


1 

• 


m , ■ 


'2n 


m 


• 


• ', 


• 


1 


2 




n 




(with L —7- oo) 






Note that all real space 





(85) 



high temperature diagrams 

up to order (tanh" /?) 
such as the closed square shown above 
are of equal value in both 
this one dimensional system 
with interaction ranges R = l,n 

and the full (L x L) two dimensional 

Ising system with nearest neighbour 

interactions 

FIG. 1. Lacing of a "two-dimensional" slab. The width of 
the slab "n" serves as an inverse temperature axis in the fol- 
lowing sense- the wider it is, the higher order in l/T that we 
may advance towards the full two dimensional model. We may 
get an analogous "three-dimensional" slab with "two imagi- 
nary time axis" if we set m — . The path will then thread 
sites before continuing upward. 

where the index b runs over the various independent loop 
momenta {k is still merely a scalar for this one dimen- 
sional problem). Unless, for a given kb integration, the 
argument of the exponent is identically zero (i.e. corre- 
sponding to a term that would be generated in the d- 
dimensional nearest neighbor problem) an "interference 
term" results. However, such a term is down by 0{A~^) 
by comparison to the "good" noninterference terms that 
occur in the d-dimensional problem. For each given as- 
signment of {c/} along the propagator lines we may per- 
form some coefficient averages over a few of the {c/} first 
and then integrate over the loop momenta and average 
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over the remaining coefficients. For "bad" interference 
terms, i.e. in those cases in which the argument of the 
exponential is not identically zero prior to the integra- 
tion, the canonical integral 



1 

2A 



dcexp[ick] 



-A 



TrSjk) 



(86) 



will result in the A ^ oo limit. When such a delta func- 
tion is integrated over the momenta in the loop integrals 
the resulting term is 0{A~^) down by comparison to a 
noninterference term for which the integral would read 



1 

2A 



dc{l) 



(87) 



-A 



In the limit A ^ oo such "bad" interference terms will 
evaporate both in the connected diagrams (for the free 
energy calculation) as well as for the disconnected dia- 
grams (included in the evaluation of the partition func- 
tion). 

Q.E.D. 

Corollary: 

For a single spin 1/2 particle with an action 



S 



dr / dr' ij{T)V{T-T')ij{T') 



= 'lp{-UJn)i'{uJ„)lp{uJn) 



(88) 



where ipiT) is the two component spinor, the averaged 
partition function is identically the same as that of the d 
dimensional nearest neighbor ferromagnet. 

The proof of this statement trivially follows from 
breaking up the [0, /3] segment on the imaginary time 
axis into L pieces and allowing L ~* oo. 

We have mapped the entire three dimensional Ising 
model onto a single spin 1/2 quantum particle! 



Employing the Hubbard Stratonovich for the quantum 
case we may make some formal one dimensional reduc- 
tions. 



Bosonization In High Dimensions By An Exact 
Reduction To One Dimension 

Fermionic electron (or other) fields may be formally 
bosonized on each individual chain. The one dimensional 
band, 



vi — —2 cos qfci, 
1=1 



(89) 



is now a simple sum of d cosines (au lieu of the standard 
single tight binding cosine). For each given set {ci}f^-^, 



the band dispersion of Eqn.(|89|) may be easily linearized 
about its two respective Fermi points. Consequently the 
standard bosonization methodology may be applied. Av- 
eraging over {ci}f^i with various weights (corresponding 
to the different observables) leads to the corresponding 
d-dimensional quantities. 

Jordan Wigner Transformation 

On mapping the three dimensional Hcisenberg model 
to a spin chain (more precisely average over spin chains) 
with sine like interactions, the Jordan- Wigner transfor- 
mation may be applied. On the chain we may set 



1 - 2a] a; 



- 11(1 - 2a]a,>l 



(90) 



where af = (erf ±iaf)/2. and the {a;} operators satisfy 
Fermi statistics 



{fli, a]} = (5y , {a^, Qj} = {a|, a]} = 0. 



(91) 



Though now the Jordan Wigner can be effected to any 
translationally invariant high dimensional problem, the 
resulting fermion problem is in general very complicated. 



VII. PERMUTATIONAL SYMMETRY 



The spherical model (or 0{n — > oo)) partition function 



Z = const H 



1 



(92) 



where the chemical potential /i satisfies 
d'^k 1 



(27r)'^^(fc) + ^' 



is invariant under permutations of {{}(fc)} 
the above, the permutations 

{h}ti - {Pk} 



(93) 

{v{Pk)}. In 
(94) 



correspond to all possible shufflings of the N wavevectors 

ki . 

This simple invariance allows all d-dimensional trans- 
lationally invariant systems to be mapped onto a 1- 
dimensional one. Let us design an effective one dimen- 
sional kernel 14// (fc) by 



5[v{k) ~ v]d'^k = I 



dVe 



eff 



dk ^y^ifik)^^- 



(95) 
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The last relation secures that the density of states and 
consequently the partition function is preserved. For the 
two-dimensional nearest- neighbor ferromagnet: 



ci 



dk 
dx 



PiVeff) 



/o VT 

and consequently 



Vl-(K//+a;-2)2 



(96) 



HVeff) = Ci 



Ve 



F(sin 



(3 — u)u ' 



-)du, (97) 



where F{t, s) is an incomplete elliptic integral of the first 
kind. Eqn. (|9^) may be inverted and Fourier transformed 
to find the effective one dimensional real space kernel 
Vef f{x). We have just mapped the two dimensional near- 
est neighbor ferromagnet onto a one dimensional system. 
In a similar fashion, within the spherical (or equivalently 
the 0{n oo)) limit all high dimensional problems may 
be mapped onto a translationally invariant one dimen- 
sional problem. It follows that the large n critical ex- 
ponent of the d dimensional nearest neighbor ferromag- 
net are the same as those of translationally invariant one 
dimensional system with longer range interactions. We 
have just shown that a two dimensional 0(n ^ 1) system 
may has the same thermodynamics as a one dimensional 
system. By permutational symmetry, such a maping may 
be performed for all systems irrespective of the dimen- 
sionality of the lattice or of the nature of the interaction 
(so long as it translationally invariant). We have just 
demonstrated that the notion of universality may apply 
only to the canonical interactions. 

The lowest order term breaking permutational sym- 
metry in our high temperature expansion is rf {x)Tf {y) . 
Thus permutational symmetry is broken to 0{(3'^) for fi- 
nite n. For a constraining term (e.g. 
for 0(1) spins) symmetric in {ri{k)} to a given or- 
der, one may re-arrange the non-constraining term 
Y.kV-^(k)\v{k)? = T,k^~HPk)\viPk)\^) and rela- 
bel the dummy integration variables H[{r](k)}] — > 
H[{r]{P^^k)}] to effect the constraining term augmented 
to a shuffled spectra v{Pk). 
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where K are reciprocal lattice vectors. 
Thus at each vertex of order 2m we have 

~5(gln-i?,0) 



(2m)! ^ 

K 



(99) 



where qin denotes the net momenta flowing into the ver- 
tex. As for all lattice models the propagators are periodic 
v{q) = v{q + K) for each of the A'' reciprocal lattice vec- 
tors K. As such the momentum conservation constraint 
may be replaced by 



1 



(2m) 



■ ^ 5(<f„, - K, 0) 



N- 



(2m; 



-(5((f,„,0). (100) 



In the aftermath, each vertex carries a factor of the vol- 
ume (A^) multiplying the discrete Kronecker delta func- 
tion. After the momentum conservation constraints are 
taken care of at each vertex, we will be left with a summa- 
tion over all remaining independent loop momenta {kh}. 
Each one of these summations is to be replaced by the 
standard 



E 



N 



Thus, at the end we will be left with a factor of 



loops -f vertices - propagator lines 



(101) 



(102) 



multiplying an integral over the independent loop mo- 
menta {kb] in Eqn. (pi]). Compounding all together, as 
in each closed diagram the number of independent loops 
-I- number of delta function momentum conserving ver- 
tices = 1 -|- number of propagator lines, a net factor of A'' 



12 



should be attributed to each independent closed bubble 
just as it does in the conventional continuum theories. 
[11] As noted earlier, for any lattice model 

v{k) ^v{k + K), (103) 

as exp[iK ■ a;] = 1 in 

v{k) = ^ V{x) exp[ik ■ x]. (104) 

The sum over x spans all discrete lattice sites. 
The way to analytically continue v{k) when ki becomes 
complex is to make the k^ dependence also invariant un- 
der shifts by Within each Brillouin Zone the k2 
dependence stems from the analytic continuation of the 
functional dependence on the real part. The only compli- 
cation is the Gibbs' phenomenon in kg which occurs when 
taking the Fourier transform of the lattice V{x). This is 
slightly avoided for the real part of V{x) as follows: Un- 
der parity and all individual reflections: V{xi,X2,X3) — 
V{xi,X2, —X3) etc. This implies that in the Fourier trans- 
form v{k): v{k\,k2,k-i) — v{ki, k2, —k^). Thus, taking 
only the even part in fca of cos(fe2 + iAfcs) we will only 
be left with the terms cos k2 cosh fcs. Now, cosh fea is con- 
tinuous and periodic when the Brillouin Zone is chosen 
as —TV < fca < TT. No Gibbs phenomenon occurs. The left 
and right derivatives are still discontinuous at the zone 
edges. 

[12] It might be easier to compute the terms in the /3 ex- 
pansion for a three dimensional Ising nearest neighbor 
ferromagnet by going to a triangular lattice and delet- 
ing all bad "constructive interference" loop integrals. A 
line threading the triangular lattice will connect a spin 
with spins of distance Range — 1, n, and n + 1 along 
the line. "Bad interference terms" can result from loop 
integrals of the form {exp[ik] exp[infc] exp[— i(n + l)fc])*'. 
This is the only way in which loop integrals can conspire 
to give nonzero results when n — > 00. Thus, trivially, the 
high temperature expansion for d = 3 cubic ferromagnet 
is equivalent to that of the triangular ferromagnet up to 
order /3^~^. Similarly, if we set s = n — 1 we obtain a 
2-Dimensional lattice model in which each spin interacts 
with its four nearest neighbors and four diagonal neigh- 
bors. 

[13] Employing f3f = -{In Z)/N = -2cosh/3 for the free en- 
ergy density in the canonical one dimensional case triv- 
ially leads to the high temperature coefficients: 



x[(ci - Xi){ci + Xi)]~ 

+5xi,o[ic2 + X2) sin7r(c2 - X2) 
+ (C2 - 2:2)] sin-7r(a;2 + C2)] 

x[(c2-a;2)(c2-f X2)]'') (106) 

- the Fourier transform of 

■C(fc) = — 2(cos cifci -I- cos 02/52), (107) 

then we will find that the average over ci and C2 of the 
coefficients in the resulting high temperature series for 
this two dimensional long ranged kernel will have the 
coefficients of the high temperature expansion of the two 
dimensional nearest neighbor Ising model. 

[14] Note the unmentionably obvious: for antiperiodic bound- 
ary conditions (as in a fermionic problem along the imag- 
inary time T axis) a Fourier space expansion is again 
possible however only the Matsubara frequencies /3uin — 
2n{n + 1/2) with integer n will appear. 

[15] The Coulomb charges take on the role of the nearest 
neighbor interactions on the original d dimensional lat- 
tice. If one would naively thread a line through the d- 
dimensional lattice (such that net number of points and 
the spin spin interactions are preserved) then one would 
find integer long range interactions on that line. An in- 
teger range interaction becomes a shifted (by the same 
range) sinusoidal Coulomb like interaction in the contin- 
uum. The last statement follows from Fourier transforms. 



This is the average of the coefficients in high temperature 
expansion in the kernel is the long ranged one of Ecm-(§) 
with d = 1. Similarly, if we set d = 2 in Eqn.(p2|) and 
average we will find the high temperature coefficients of 
the two dimensional nearest neighbor Ising model. 
Just as trivially, if we examine long ranged interactions 
in two dimensional spawned by 

V{xi,X2) = -2^5i,2,o[(ci + Xi) sin7r(ci - xi) 
+ {ci — xi)j sin7r(a;i -I- ci)j 
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